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Quantum mechanical averaging of the particle concentration operator is an effective starting point 
for derivation of the many-particle quantum hydrodynamic equations. In many-particle quantum 
systems, we have to separate the ordered motion of the local center of mass (velocity field), the 
thermal, and the quantum motion. The quantum mechanical average process, invoked here, is com- 
pletely determined, and is different from the usual averaging processes that introduces undefined 
probabilities for quantum states. It is shown that the Madelung decomposition for the iV-particle 
spinor wave function allows the correct introduction of the velocity field, and gives explicit expres- 
sions for the quantum contributions to both the momentum, and the spin flux. The formalism 
also contains plasma effects produced by the Coulomb and spin-spin interparticle interactions. It is 
shown that both interactions appear in plasma dynamics as effective electric and magnetic fields. 
As result we find a fully coherent description for spin quantum plasma in the self-consistent field 
approximation for interparticle interactions. A simple consequence - the change brought about by 
interparticle correlations on the propagation of electrostatic waves in a spin quantum plasma - is 
discussed. 

PACS numbers: 52.25.-b, 67.10.-j, 67.30.hj 
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I. INTRODUCTION 

When a one-particle Schrodingcr equation is expressed 
as a set of two "hydrodynamical equations" : the continu- 
ity equation and the Euler equation, [l|, , it differs from 
the classical Euler equation in that the thermal pressure 
in Euler equation is replaced by an explicit quantum pres- 
sure whose origin is the quantum Bohm potential. The 
latter, proportional to the Ti , is caused by the nonlocal 
nature of the particle formed by a de Broglie wave packet. 

For scalar particles, only two real classical functions 
(an amplitude defining the one-particle quantum fluid 
density, and a phase that characterizes the one-particle 
velocity as a potential field) are needed to represent 
the wave function. This decomposition produces the 
quantum Bohm potential [ll, 0- For spin-1/2 particles, 
obeying Pauli-Schrodinger (P-S) equation, however, the 
spinor wave functions require a generalization [3J, |j| , that 
give rise to new spin contributions to quantum pressure, 
as well as a spin evolution equation. 

While the "fiuidization" of one-particle Schrodingcr 
and P-S systems is straightforward, the quantum ensem- 
ble average that will lead to a desirable macroscopic fluid 
description for a spin quantum plasma has to be carefully 
performed. In the often quoted works (see for example 
Refs. Q, Q); the ensemble average is understood to be 
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over the normalized occupation probabilities for a given 
mixture of quantum states. In this form, each macro- 
scopic quantity is defined as the sum over the correspond- 
ing fluid-like elements modulated by its probability of oc- 
currence. The macroscopic equations, so obtained, have 
been investigated for a host of new interesting linear and 
nonlinear quantum-magnetohydrodynamical phenomena 
(see, for example, Refs. |5l4l7| and references therein) 
including new instabilities driven by spin [17H22| . A 
new revealing formulation casts the system into a vor- 
tical structure where the the spin-dynamics is shown to 
resemble the vortical-dynamics of an ideal fluid [23[ ; ex- 
tensions to relativistic quantum plasmas have also been 
constructed [24J| . 

Although the procedure developed in Refs. 0, Q may 
be consistent for non-interacting particles, it is difficult 
to see how it could deal, properly, with 7V-wavcfunctions 
correlations (as those produced by Pauli exclusion prin- 
ciple), and inter-particle interactions. In this work , we 
seek to develop a more consistent quantum mechanical 
averaging process to derive many-particle quantum hy- 
drodynamic equations. The system is defined through a 
many-particle Schrodingcr equation with a many-particle 
wavefunction ip{r\, ..., rjy, t) [shortly i/j(R, t)], where r p is 
the coordinate where the p-th particle is located; R des- 
ignates the whole coordinates set. Thus, we are able to 
introduce Coulomb and spin-spin interparticle interac- 
tions constructing a fully coherent description for spin 
quantum plasma in the self-consistent field approxima- 
tion. One of the tests that the new formalism must pass is 
that it will produce the corresponding one-particle limit 
in the approximation of independent particles. Similar 



attempts in kinetic theories for quantum plasmas has 
been proposed [25j . 

The paper is organized as follows. In Section [ill we 
present the basis for the macroscopic equations using the 
quantum averaging for the iV-particle Pauli equation. In 
Sec. HIT] we obtain the fluid equations using the general- 
ized Madclung decomposition for the wave function. The 
fluid version of the iV-particle Pauli equation will contain 
the contribution of the spin, producing additional forces 
in the Euler equation and an evolution equation for the 
magnetic moment. Besides, the momentum flux and the 
spin flux of both equations contain quantum and thermal 
effects. We show that the inter-particle forces produce ef- 
fective electric and magnetic fields that interact with the 
fluid velocity, and with the spin fields. In Sec. II VI we dis- 
cuss the approximations needed to obtain a close form for 
the Bohm potential, and in Sec. [V]we give the vortical 
description for the present formulation. In Sec. IVII we 
study the contribution of the previous mentioned effects 
for the propagation of electrostatic waves. And in Sec. 
IVIII we discuss our results. 



II. CONSTRUCTION OF MACROSCOPIC 
EQUATIONS 

The dynamics of the iV-intcracting fcrmions is given 
by the Pauli-Schrodinger equation 



ihd t ^(R,t) = H^(R,t), 



(1) 



with the Hamiltonian 



H = 



N 

E 



D l D l 

2m pp 



fia l p B; 



E 

p,n,p<n 



H-o 



(2) 

acting on the ./V-particle wave function ^(R, t). Ein- 
stein convention over repeated latin indices will be used 
throughout this paper. Here m and q are, respectively, 
the mass and charge of the particle, \i = qfr/(2mc) is the 
magnetic moment, h is the reduced Planck constant and 
c is the speed of light. The covariant derivative operator 



Dl 



-thai, 



-A 1 
, A p- 



(3) 



and cj> and A 1 are the scalar and vector potential, o % are 
the Pauli matrices and B % are the components of the 
magnetic field. The last term in the Hamiltonian @ 
describes the long-range interaction between particles in 
a self-consistent field approximation [4|, [28[ 



H pn = q G pn — fi G^ n a p a, 



3 
n > 



(4) 



where the first term represents the Coulomb interaction 
and the second term is the spin-spin interaction. The 
two Green functions arc: 



(-Jpn — ***pn\\^p ^-n\) 



1 



(5) 



and 0, [H 



Gil = Gl 



n\\ x P 



= d % 8> 



r»l) 

1 



Tr, Yf] 



+ 4n6 ij 6(r p - r„) . (6) 



The quantum mechanical average that will yield the 
macroscopic variables, can be defined in a standard way 
usingthe iV-particle wavefunction. The microscopic den- 
sity |4l. l2q. I29m2| . for instance, will take the form, 

N 

n = n(r,t) = / dR^, 6(r-r p )tf(R, t)^(R, t) = (* + #) , 
J P =i 

(7) 
with dR = Ii p=1 dv p . In the rest of the paper, (£) = 

/ dR J2 p= i ^( r — r p)C wu l symbolize the average of the 
quantity £(R,t). 

By taking the time derivative of (O, we derive the 
first conservation law, the continuity equation: 



£ + *' = <>■ 



(8) 



where the macroscopic current (h.c. stands for hermitian 
conjugate) 



J< = J l (r, t) = ^- (&Di* + h.c.) 
2m F ' 



(9) 



is the appropriately averaged microscopic current. Fol- 
lowing Ref. [J], the evolution equation for the current 
,P is obtained by taking the time derivative of © and 
invoking the P-S equations (fTJ) and ([2]) 



dJ l 
dt 



-9,n ,J = ^-nE i + -^- £ ljk J j B k + — M k d l B k 
m m mc m 

_«_ fdr'd l G(r,r')n 2 (r,r\t) 
m J 

+ 1 /'dr'G)' ; G''' fc (r,r')A/f(r,r',i) I 
m J 

(10) 

where E l = E l (r,t) and B l = B l (r,t) are the macro- 
scopic electric and magnetic fields. Four macroscopic 
density tensors are introduced in (fTU|) . The first one is 
equivalent to the momentum flux in fluid dynamics, 

IF = U ij (r, t) = ^- (¥D p D p V + (Z?;*) f Djtf + h.c.) 

(11) 

and the second one is the magnetic moment density 

M l = M*(r,t) =/x(*V**) . (12) 

The other two new fields are related to inter-particle 
interactions. The two-particle number density (corela- 
tion function) 

II 2 



(r,r',f)= j dR J2 S{r-r p )5(r'-r n )^(R,t)^(R,t), 

(13) 



p,n,p<n 



is created by Coulomb correlations, and the tensor 
Mi j (v,r',t) = 



(14) 



p,n,p<n 



and spin A lJ fluxes. We will soon see that the Madclung 
decomposition (j2"0"j) does, indeed, yield the correct ex- 
plicit form for IT- 7 and A y . 

We first, using Eq. ([7]), see that the macroscopic den- 
sity 



<« 2 >- 



(21) 



is the manifestation of the spin-spin interaction between 
two particles. 

In order to formally complete the particle dynamics, 
we must also calculate the dynamical equation that gov- 
erns the evolution of the macroscopic spin magnetization. 
Taking the time derivative of the magnetic moment den- 
sity, we derive [4| 



dAP 



■W 



^t e ijk M j B k 



2p 



'jk 



: f dr'G km (r,r')M^ m (r,r',t), 

(15) 
where the new magnetic moment flux tensor stands for 

A ij = A l] (r,t) = #- (^aiDiV + h.c.) . (16) 

2m F F 

Finally, the whole system must be coupled to the 
Maxwell equations 



diB* = , 

e l3k diE k 



d % E l = iirqn , 
= --d t B\ 



e ijkg jgk = ^Q_ji + i7re ijkg M k + lg E i 

c c 



(17) 
(18) 

(19) 



where the total current is formed by the charged fluid 
current J 1 and the magnetization current e % ^ k djM . 



III. MACROSCOPIC FLUID DESCRIPTION 

We now proceed to construct a hydrodynamical theory 
from the macroscopic equations derived in the last sec- 
tion. The first step is to do a Madelung decomposition of 
the TV-particle wave function (see, for instance, Ref. [371 ]) 



(20) 



*(i?, t) = a(R, t) exp ( tS(R, t) ) p{R, t) , 



where a is an amplitude, S is a phase and tp is a normal- 
ized spinor (ip'ip = 1). 

We expect that the Madelung decomposition, first in- 
voked for spinless particles, will help to obtain an explicit 
form of the quantum Bohm potential 
out the decomposition used in Refs. 
particle wavefunction of spinning particles failed to pro- 
vide the correct convective part of the momentum IP- 7 , 



26IJ27J. We point 
23 for the N- 



The Madelung decomposition implies, for the p-th par- 
ticle, a microscopic velocity 



vl(R,t) 



1 



dtS 



ih 

— p'OpP 

m H mc 



' f ai 



K , (22) 



that will help write down the macroscopic current ^ as 
J 1 = (a 2 v p ) , (23) 

and a macroscopic fluid velocity as 



V l = V l {r,t) 



(24) 



n(r,t) 

In the same way, from the microscopic spin 

4(fl,t) = ¥>V*¥>, (25) 

normalized to unity, we derive the magnetic moment den- 
sity ([12]) 



M*= M (a 2 4>, (26) 

such that the macroscopic spin density emerges as 

M l (r,t) 



S l (r,t) 



Hn(r,t) 



(27) 



Let us now evaluate the momentum flux (jlip in terms 
of the newly defined fluid variables, 



(28) 



^ = {^( d > d l a - ad P d » 



2m 



h 2 a 2 
+a 2 mv l p v p + —dls k di p s k 



The first two terms of the tensor create the Bohm po- 
tential, the third contributes to the fluid pressure, and 
the last term is a force produced by the self-interaction 
of the spin particle. The magnetic moment flux density 
(JTHJ) is, similarly, expressed as 



= n(a 2 8* p v> p 



2m' 



Jmfc „m 



s m d J s k 
p p p 



(29) 



As a further step in the "hydrodynamic" project, we will 
introduce the normal as well as the spin temperatures. 
Let us split 



z l = v l - V 1 



"'n 



S\ 



(30) 



such that {a 2 z p ) = and {a 2 w p ) — 0. The field z p signi- 
fies thermal speed, and w p is the displacement (thermal) 



fluctuation of the spin about the macroscopic average 
S l . Both displacements will be used in the fluxes ([25)1 
and (|29| . Note that the macroscopic spin is no longer 
normalized to unity because of the thermal-spin effects, 



O Oi 



1 - -(a 2 w' p w p ) 



(31) 



Because the evolution of the fluid equations is dom- 
inantly driven by the electromagnetic field interaction, 
we could perform a self-consistent approximation for the 
correlations (jl 



We are ready now to write down the dynamical equa- 
tions for the spin quantum plasma explicitly in term of 
the fluid variables: we begin with the continuity 



On 
It 



+ d t (nV l ) = . 



(33) 



ra 2 (r,r',i) « n(r,t)n(r',t) , 
Mi j (r,r',t) « M i (T,t)M j (r f ,t). 



(32) the momentum evolution equations 



a + ™,iv 



—E l + -L e ijk V j B k - -dip ij + ^S j d i B j + -^<9 4 (d i S k d i S k ) - —Q l , 

8m z 



in 



mc 



m 



nit) 
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(34) 



where we have used the relation (|A1|) . We have also 
introduced new definitions: the pressure tensor p lJ = 



(a 2 z l p z 3 p ), the effective electric 



E i = E l -qd l / dr'G(r', r)n(r', t) , 



(35) 



fluid pressure originating in the self-interaction of spin 
fluctuations. 

Finally, the spin evolution equation takes the form 



d_ 
di 



V 3 d 3 ) S l 



^t e ijk S j B k --Q i , (38) 

n n 



and the effective magnetic field 

~. . he f 

B l = B l + dj (n&S^+fi / dr'G y (r , ,r)n(r',i)5 J (r',i) 

2qn J 

(36) 

The last term in (J34]) is the quantum force produced by 

the Bohm potential and the thermal-spin interactions 



where 



^ = M rf K fl ) - ^ d > + -^^4 



2m 
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—dlstdlwt - —dl.w'idlw 
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di 



4 d >p 



m 
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in 



„ k d J w k 

p ^p^p ^p 



h 2 n 
8m' 



d i 



; d i 



& (1 <a a ™*^>) 



(37) 



and it is rather complicated. 

In Eq. (|3"4"|) . the effective electric field ([33)) contains the 
intrinsic electric field generated by Coulomb interaction, 
which is, naturally, the gradient of a potential. The effec- 
tive magnetic field ([55)1 is also enhanced by the magnetic 
field generated by the variations of the spin field, and 
by the intrinsic magnetic field produced by interparticle 
spin-spin magnetization. In this sense, the present ef- 
fective magnetic field is a quantum generalization of the 
one used in Ref. [2J]. The fifth term on the right-hand 
side (rhs) is a perfect gradient and acts as an additional 



Q l = d. 



a 2 z j s i - —£ imk s m d j w k 
p p m p u p w p 



(39) 



stands for all the thermal-spin interactions. To derive 
([55]). we used (|X2"j) . From Eq. (J3HJ) we note that the 
macroscopic spin is normalized only if the quantum spin- 
thermal interaction Q is neglected. This is consistent 
with (|3Tj) . 



APPROXIMATE FORM OF QUANTUM 
BOHM POTENTIAL 



IV. 



Due to many-particle effects, the force Q l (|37| appear- 
ing in Eq. (|3"4"]) j26l . |27| is rather complicated, and con- 
tains implicit averages that defy an expression in terms of 
hydrodynamic variables, the particle concentration, the 
velocity field and the average spin field. In order to feel 
confident that the quantum pressure, though, somewhat 
opaque, is really correct, we would show that in the sim- 
plest approximation (noninteracting particles) , it reduces 
to the correct known limit of the quantum pressure in 
one-particle QHD, the Bohm potential. 

We will concentrate on the first two terms of Q l , which 
can be extracted out of the averaging operator in the ap- 



proximation of weakly interacting particles. This trans- 
forms the quantum force in (|37l) to 



h n 



2m 



A x 



+-8, (a 2 {di p w k p d;4 + 44 9 >p " d> k P d> k v )) 



h n 



V -(»%') 



(40) 



that explicitly contains the correct Bohm potential in the 
first term on the rhs. 



V. VORTICAL STRUCTURE 

If thermal-spin interactions are neglected, the quan- 
tum force (|4U|) becomes a perfect gradient of the Bohm 
potential, then the system of equations ([55]). (|M|) and 
(f38|) are reducible to a vortical structure that looks ex- 
actly like the one developed in J23[. The principal re- 
sult of Ref . [23( , dwelling on the standard formulation of 
the Pauli-Schrodingcr equation, is the construction of a 
classical-quantum hybrid vorticity 



m = m 



2q q 



(41) 



from the classical £V C = B l + (mc/q)e' i: ' k djVk, and the 
quantum part fi* = (l/2)e ijk e lmn Sid j S m d k S n [H- The 
hybrid vorticity obeys the Helmholtz vortical dynamics 
leading to a new conserved hybrid-hclicity 



h_ = 



d 3 x Q* 



P'_ 



(42) 



such that dh-/dt = 0. Here, P l _ is the generalized po- 
tential defined by £l % _ = s l: > k djP k . Notice that B l given 
in (|36[) contains the extra term (the Green function in- 
tegral representing correlations) , and thus generalizes all 
previous results including those of Ref. [23j . For more 
details on the spin quantum vortex dynamics, we refer 
the reader to 1231] - 



VI. ELECTROSTATIC WAVES 

The intcrparticlc Coloumb and spin-spin interactions, 
cause, inter alia, interesting modifications to classical 
wave propagation. As an illustration, we will calculate 
the electrostatic waves in a two-dimensional plasma with 
isotropic pressure p. In this case, the Green function ^ 
has been shown to be [3a| 



G ij = 

where I = \l\ 
component of I. 



S ij - WV jl 2 8tt 



S ij S(l) . 



(43) 



P 3 

- r„|, and l l stands for the i 



Let us explore a simple model in which the two- 
dimensional plasma is composed of electrons {q = — e; 
/i = —eh/2mc < 0) moving in a plane perpendicular 
to the external constant magnetic field B$e z and with 
pV = S^p/m. To study electrostatic waves, the total 
fluid density fluid is taken to be n = no + 5n, where no is 
the background density and 6n <C no is the perturbation. 
In the absence of an equilibrium flow, the perturbed fluid 
velocity, 6V 1 , is such that 5V Z = 0. The perturbations 
of the macroscopic spin field 5S l are perpendicular to 
the normalized equilibrium spin field, that is aligned an- 
tiparallel to the external magnetic field to minimize the 
magnetic moment energy. All perturbed quantities vary 
as explikjT^ — iuit), where r* and k % are restricted to the 
x — y plane. The perturbed continuity equation simplifies 
to 



Sn 



noki 8V l 



(44) 



while (neglecting the thermal-spin interactions) the lin- 
earized momentum equation yields 
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(45) 



= eBo/mc is the cy- 
is the thermal 



where SE l is the electric field, 
clotron frequency and v s = y/ksT/m 
speed. For simplicity, we consider electric field perturba- 
tions through Poisson equation (however to include the 
intrinsic electric field in ([33)) is straightforward). The last 
term, representing the effect of correlations, contains the 
function 



/3{k) = 2tt 



kl„ 



dx 

— Jo{x) 



(46) 



where Jo is the Bcsscl function, of order and l m in is 
the minimum distance for the spin-spin interaction. This 
integral appears from the spin-spin interaction in the def- 
inition (|3"6"| for the effective magnetic field (see Appendix 
[Bjl . It must be evaluated numerically. 

In addition, the perturbed spin evolution equation be- 
comes 



iioSS 1 



V Ti n 2m 
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2tt 



dl 2n 



(47) 



, (48) 



T jk = \ dr'G jk (r',r)e 



ik-fr 7 — r) 



(49) 



which is evaluated in (|B4|) . 

Finally, the perturbed Poisson equation for a two- 
dimensional plasma reads [36| 



ikj 8E l = —2nekdn . 



A. Spin waves 



(50) 



Although the exact value of /3(k) requires numerical eval- 
uation, approximate values may be trivially calculated 
for small as well as large kl m i n . Restricting to small 
kl m in, we find /3(k) ~ 2TrJo(kl mm )/kl mm ~ 2?r /kl min . 
The correlation term, then, goes as k 2 and competes with 
the thermal term in providing quadratic dispersion. Simi- 
lar results have been found for plasmas with dipolc-dipolc 
interactions 1271. 



Because the modes are electrostatic, the spin equations 
(|47|) are decoupled from the rest. Consequently there 
are two independent branches of the dispersion relation 
that can be labelled, respectively, the spin waves and the 
plasma waves. The spin wave dispersion relation 



u 2 = 



4irfi 2 nok 



where 



Q 



2[iB hk 2 An^no 
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(51) 



(52) 



can be further simplified in the "long" wavelength limit 
2-Kklmin <C 1, i.e, when Airy^nQk/fi <C Q; being the 
length l m in ~ n-o -1 ' 2 , the interparticle distance in this 
two-dimensional model. The resulting dispersion relation 



hk 2 4iTfj, 2 no 

ZTH lilm.in 



(53) 



is a generalization for the dispersion of spin waves found 
in Refs. [38|, |39|. The first term on the rhs of (|53|) is 
the cyclotron frequency contributed by the external field, 
and the second one is the standard quantum contribu- 
tion. It is the third term, originating in the magnetic 
field induced by the spin self- interaction (|36[) , that is new. 
The new quantum correction (coming from multi-particle 
correlations) becomes comparable to the standard term 



When (klmin) ~ ^classical/ h 

classical electron radius. 



where A 



do 



is the 



Plasma waves 



The dispersion for the "plasma branch" (in the two 
dimensional plasma) can be found using Eqs. (|44|) . (|45] 
and (|50j) . 
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(54) 



In this hybrid mode, the first term is the effective plasma 
frequency in a two-dimensional plasma |36| . the second 
term is the cyclotron contribution, and the third one 
comes from the thermal effects. The two last terms are 
the quantum contributions: the fourth term is the well 
known contribution of the Bohm potential, while the fifth 
one, proportional to /3 = /3(fc), is a totally term originat- 
ing in contributions from the spin-spin interaction (1461) . 



VII. CONCLUSIONS 

Starting from the A-particle Pauli-Scrodinger equa- 
tion, we have derived the fluid equations for a spin quan- 
tum plasma, using a well-defined quantum mechanical 
averaging procedure. This rather general averaging pro- 
cedure, in addition to providing proper definitions for 
macroscopic fluid variables, allows us to introduce inter 
particle corclations caused by Coulomb as well as spin- 
spin interactions. 

We would like to stress that the ability to express the 
interparticle interactions in terms of fluid variables is a 
fundamental strength of this class of quantum averag- 
ing procedures. Consequently, in addition to the usual 
quantum forces due to the Bohm potential and the par- 
ticle spin, new forces due to interparticle interactions ap- 
pear; the latter can be viewed as effective electric and 
effective magnetic fields. The structural similarity of the 
new forces with the electromagnetic ones, makes it pos- 
sible to cast the whole quantum plasma dynamics into 
a generalized vortex dynamics system with a conserved 
generalized helicity that has additional content owing to 
interparticle corelations. 

As an application of the formalism, we calculated the 
effects of interparticle interactions on the propagation 
of electrostatic waves; the dispersion characteristics of 
both the spin and the plasma-cyclotron mode can be pro- 
foundly changed. 

Finally, we want to emphasize that the operational do- 
main of the very straightforward quantum averaging pro- 
cedure presented in this work, can be readily extended to 
deal with various possible particle interactions that can 
take place in a quantum plasma. 



Appendix A: Useful relations 

The following relations are useful in previous deriva- 
tions 

dj (nd :i S k d l S k ) + ^d l (d j S k d j S k ) 

+nS k d l (-dj (n&>S k )) = nd i (-d k (nS : >d k S : >)\ , 



(Al) 



s l]k S,d m (nd m S k ) = d m (e^nSjdmSk) ■ (A2) 



Appendix B: Green functions 

The following integral 



dr'G(r',r)e ik < r '-^ = d$ 






(Bl) 



with £ = r' — r and £ = |£|, is readily evaluated for 
particles constrained to a two-dimensional plane, 
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f dr'G{r' 7 r)e lk(r '- r) ■- 
Using the first part of (|43|) . the integral 
j dr'Gy(r',r) e jk ( r '- r ) = / ^ G„ n (4)e ik « , (B3) 



is simply 
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ik« 



(B4) 
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